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Abstract 

A two-parameters family of Backlund transformations for tlie classical elliptic 
Gaudin model is constructed. The maps are explicit, symplectic, preserve the 
same integrals as for the continuous flows and are a time discretization of each 
of these flows. The transformations can map real variables into real variables, 
sending physical solutions of the equations of motion into physical solutions. 
The starting point of the analysis is the integrability structure of the model. It 
is shown how the analogue transformations for the rational and trigonometric 
Gaudin model are a limiting case of this one. An application to a particular case 
of the Clebsch system is given. 
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1 Introduction 



The Gaudin models describe completely integrable long range spin-spin systems, both 
at the classical and at the quantum level. They were first introduced by Gaudin [7], as 
the anisotropic, integrable generalization of what today is called the xxx, or isotropic, 
Gaudin model; the results of the analysis were the construction of the xxz model, for 
the partially anisotropic case, and of the xyz model, for the fully anisotropic case. The 
Lax matrices depend on the spectral parameter respectively through rational, trigono- 
metric and elliptic functions. The Poisson structure of the models can be specified in 
the framework of the r-matrix approach [22]; in [21] Sklyanin showed how to obtain 
both the Lax matrices and the r-matrices of the models by a limiting procedure on the 
lattice Heisenberg magnet. Obviously the integrability structures of the trigonometric 
and rational cases can be also obtained as particular realizations of the corresponding 
structures of the more general elliptic model. From the point of view of separation of 
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variables, functional Bethe ansatz and quantum inverse scattering method, a lots of 
results appeared in the late 80s-early 90s [10] [9] [23] [28] [29]. On the other hand a 
great number of works have pointed up until now very interesting connections between 
the Gaudin models and various branches of physics. For a recent list of the subjects 
involved the reader can see for example [TSj. Here, in order to give an idea of the 
heterogeneity of the concerned issues, we wish only to mention the connections with 
the BCS theory and small metallic grains [1], with pairing models in nuclear physics 
[6], with the Coulomb three-body problem [12], with Lagrange and Kirchhoff tops [T6] . 
In [8], [25], [26] the authors gave the Backlund transformations for the rational Gaudin 
model in the framework of a larger research program started some years before by 
Kuznetsov and Sklyanin (see [H]) on the applications and properties of such transfor- 
mation to finite dimensional integrable systems. In particular in [8] it was shown how 
the two-parameters family of transformations obtained can be seen as a time discretiza- 
tions of a family of flows and that the interpolating Hamiltonian flow encompasses all 
the commuting flows of the model. The features of the transformations found here are 
noteworthy: in fact the maps are (i) explicit, (ii) symplectic, (iii) they preserve all the 
integrals of the continuous model, (iv) they possess the so called spectrality property 
[13] |14] . Furthermore they commute, that is the composition of two different trans- 
formations does not depend on the order of application, and with a two-parameters 
transformation it is possible to send real variable into real variables (so that the Back- 
lund transformations send a real solution of the equations of motion into another real 
solution). 

In the wake of these results in 2010 the analogue constructions for the trigonometric 
case were found [19] [20] . All the above properties hold true in this case too. In [20] the 
authors have also shown how, in the limit of small angles, the transformations obtained 
give exactly the results of the rational case as given in [S]. 

The aim of this work is to complete the picture finding the Backlund transformations 
and the corresponding time-discretization for the elliptic Gaudin model. The paper is 
organized as follows: in section ([2]) the main features and the integrability structures 
of the elliptic Gaudin model are briefly revised, in section ([3]) the dressing matrix is 
constructed; the parametrization obtained gives directly the dressing matrix of the 
trigonometric case [20] when the elliptic modulus k of the Jacobi elliptic function is 
set equal to zero. In section (j4]) the explicit form of the transformations is found by 
using the spectrality property. The symplecticity of the maps is also discussed. For any 
fixed set of initial conditions the transformations turn out to be birational maps. In 
section it is shown how it is possible to obtain a two-parameters family of physical 
transformations, that is transformations sending real variables into real variables. In 
section iQ the continuous limit of the discrete dynamics is found: the interpolating 
Hamiltonian flow again encompasses all the commuting continuous flows of the model. 
Lastly, as an application of the results found, we construct the Backlund transfor- 
mations for a particular case of the Clebsch system [1] through a procedure of pole 
coalescence on the Lax matrix of the Gaudin model and taking advantage of the fact 
that this procedure preserves the r-matrix structure. The results are a generalization 
of the Backlund transformations for the Kirchhoff top as given in [2T] . 
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2 The elliptic Gaudin magnet 

The elliptic Gaudin model is defined by the following Lax matrix: 

m = 



>^cn(A-^ 



AiX) B{X) 
C(A) -AiX) 



^^^^^^s^-is^dniX-X, 



sn(A - Xj] 



(2) 



sn(A — Xj 



In (IT]) and ([2]) A G C is the spectral parameter, Xj are a set of arbitrary real 
parameters of the model, while (sj, sj, j = 1, . . . , iV, are the dynamical variables 
(spins) of the system. In the r-matrix formalism the Poisson structure is fixed by the 
following equivalence [28] : 



{L(A),L(/x)} = [re(A-/i),L(A)® 11+11 ®L(/i)], 



(3) 



where rg is the elliptic solution to the classical Yang-Baxter equation of Belavin and 
Drinfel'd [5] , [3] , [28] : 



sn(A) 



f cn(A) 



1 l~dii(A) 
\ 2 






l+dn(A) 
2 







l+dn(A) 
2 






l-dn(A) \ 
2 ' 




cn(A) J 



(4) 



The r-matrix structure (l3|) entails the following Poisson brackets for the functions 
{A(A),A(/i)} = 0, 



dn(A - 

sn(A — /i) 
1 — dn(A — /i) 
sn(A — /i) 

^ C(A)(1 - dn(A - /x)) - B{X){1 + dn(A - fi)) + 2S(/i)cn(A - fi) (5) 

2sn(A — /i) 

^ ^(A)(dn(A - At) - 1) + C(A)(1 + dn(A - fi)) - 2C{fi)cn{X - /x) 

2sn(A — /i) 
(A(^)-A(A))(l + dn(A-/i)) ^ 
sin(A — //) 

Equivalently, the spin variables (sj , sj, s|), j = l..A^, have to obey to the corresponding 
algebra: 



{i?(A),i?(/i)} = i(A(A) + A(/i)) 
{C(A),C(/i)} = i(A(A) + A(/.)) 
{A(A),B(/i)} = i 
{A(A),C(/i)} = i 
{i?(A),C(^)} = i 



{sj, si} = 6jksl, {sj, 4} = 6jksl, {s'j, si} = 5jkS 



(6) 
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Due to the direct sum structure of the Poisson bracket the square length of each 
spin is a Casimir function for the elhptic Gaudin model, so we have N Casimirs, given 
by 

{s^r + {s^r + is^.r = j = i..n 

Furthermore the model possesses N integrals of motion in involutions w.r.t. the Pois- 
son brackets (jH]): the determinant of the Lax matrix is a generating function of such 
integrals (see 

- det(L) = A\X) + S(A)C(A) = V + 2wH,aw{\ - A,)) - (7) 

where C is the WeierstraB zeta function, w = {ei — 62)"^, = p(^) and {wi,W2) are 
the periods of the WeierstraB p function (see|X]). The Hamiltonians Hi are given 
by: 

^ slslcn{\i - Afc) + g^g^dn(A^^ - A^) + gf 

Note that only A^— 1 among these Hamiltonians are independent, because of Hi = 0. 
The other integral Hq is given by the formula (see 

N 

Ho = Y,istsldn{Xi - \k) + k\\slcYi{X, - A,)) + 



i,k 
N 



{\i - Afc) (gf g^cn(A^ - Afc) + 5^dn(A^ - A^) + gg) (9) 



a(A) = ro(C(tx7A) - C(2wA)) ^ 



sn(2A) 

Due to the existence of an r-matrix, the Hamiltonians Hi are in involution for the 
Poisson bracket ([6]): 

{H,,H,} = Q z,j = 0,...,iV-l (10) 
The corresponding Hamiltonian flows are then given by: 

ds^- „ dsf , 

it-^"""^ it-^"-^n (11) 



3 The dressing matrix 



By an Hamiltonian point of view, Backlund Transformations for finite-dimensional 
integrable systems are (families of) symplectic maps preserving the integrals of motion 
\1'S\ : so we are searching for an expression relating the dynamical variables (■^J, gj, gj) 
to the new set of variables (g|, gj, g|) and such that the brackets and the integrals 
(IH]) are preserved. The generating function of the integrals is given by the determinant 
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of the Lax matrix ([T]), so the dressed Lax matrix L{X), that is the Lax matrix of the 
tilded dynamical variables, has to have the same determinant of L{X). This means that 
the two Lax matrices are related by a similarity transformation by means of a dressing 
matrix -D(A): 

L(A)D(A-Ao) = I^(A-Ao)L(A) (12) 

Here Aq will be one of the two Backlund parameters: actually we will construct a 
parametric family of transformations and this will be a crucial point when we will 
request to the maps to be "physical", that is to send real variables in real variables. 
The rational and trigonometric Gaudin models are limiting case of the elliptic one. 
The dressing matrices for the Lax matrices of these models are given, respectively, by 
the elementary Lax matrix of the xxx and xxz Heisenberg spin chain on the lattice [S] , 
[20] . It is obvious, for the dressing matrix of the elliptic Gaudin model, to make the 
ansatz of the elementary Lax matrix of the xyz Heisenberg spin chain. Note however 
that this matrix has to enjoy the same symmetry properties with the Lax matrix ([1]). 
In fact, the quasi-periodicity dHS]) of the Jacobi elliptic functions entails the following 
formulae for L(A) (see also [29]): 

L(A + 2K) = asL{X)a3 L{X + 2iK') = aiL{X)ai (13) 

where with cTj, i = 1, 2, 3 we indicate the Pauli matrices, and K and K' are respectively 
the complete elliptic integral of the first kind and the complementary integral ( 167|) . 
The above points suggest to make the following ansatz for D{X) 

D{X) = ^ol + + dn(A)^2a2 + cn(A)^3fX3) (14) 

sn(^A j 

This is exactly the one-site Lax matrix for the xyz Heisenberg spin chain on the lattice 
[5]. The symmetries (fT3|) are also preserved: 

D{X + 2K) =(rsD{X)(rs D{X + 2iK') = aiD{X)ai (15) 

So far, J^i, i = 0..3 are four undetermined variables, but we are free to fix one of them 
because of the homogeneity of the equation fll2p . The Lax matrix ([T]) has simple poles 
at the points A = Xj (mod 2K,2iK'), j = 1..N; the relation f ll2p is an equivalence 
between meromorphic functions (the elements of the matrices), so that we have to 
equate the residue at the poles on both sides. Because of the symmetries (|T3l) and f|T5|) 
we can look only at the poles in A = Aj, j = 1..N, that is: 



LjD, = D,L, (16) 



where 



Lj = ^^^y ^ J , D, = D{X = A,) (17) 

In principle equation (ITB]) gives an implicit relationship between the old (untilded) 
variables and the new (tilded) ones. It is however possible to get an explicit relationship 
by recurring to the so-called spectrality property [13] [H]. To this aim, one need to 
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force the determinant of the Darboux matrix D{X) to have two nondynamical zeroes for 
two arbitrary value of the spectral parameter A, say for A = Aq ± /i. This leaves us with 
only two undetermined variables in (fT^ . As we will see, the spectrality property will 
fixes these two variables, that we will call P and Q, as functions of only the untilded 
dynamical variables, so that the maps defined by (1161) will be explicit. 
Summarizing, by fixing for simplicity = 1 and imposing the constraints 



det{D{X - Ao)) 



= 

X=Xo±fi 



we are left with two undetermined parameters, that we denote with P and Q: by 
choosing a particular parametrization of the constraints we can write: 



-1- -r I^^SsniX) sn( 

sn(A) ^ *"^3sn(A) 



^(A) = I , :tUx) 1 with { 



^ cn(/i) 
i^i + ^2 dn(^) = P 
[ i^i - ^2 dn(/i) = Q{2 sn(^) - PQ) 

(18) 

We recall again that P and Q are undetermined dynamical variable and that Aq and /x 
are constants: they are parameters for the Backlund transformations. Note also that 
with this parametrization, in the limit A; — > 0, one obtains the dressing matrix for the 
trigonometric Gaudin model [19] (up to a trivial multiplicative factor inessential for 
the form of the Backlund transformation, as explained before). 



4 Backlund transformations 

Now we make use of the spectrality property to find P and Q in terms of only one set of 
variables, the untilded ones. The matrices {D{\ — Aq)) and (-D(A — Aq)) 

A=Ao+/i A=A()— /I 

are of rank one. We call and their respective kernels. By acting with these 
kernels on the equivalence defining the Backlund transformations, we see that they are 
also the eigenvectors of L{\q + y) and L{Xq — y): 

L{\^±y)D{±y)\n±) = = D{±y) [L(Ao ± /i)|fi±)] =^ L(Ao±//)|fi±) = 7±l^^±) (19) 

By viewing the generating function of the integrals ([7]) as a function of A, we define: 

^^[X) = -det{L{X)) = A{X)^ + B{X)C{\) (20) 

where ^(A), B{X) and C(A) are given by ([2]). The two eigenvalues are then given by 
7± = 7(A) . The two kernels are expressions of the variables P and Q, so 

A=Ao±/i 

that the eigenvectors relations ( IT9|) for L(Ao ± y) link these variables with the elements 
of the Lax matrix of the untilded variables. Explicitly, the two kernels are given by: 
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and these expressions in turns lead to the formulae: 



B(A) 



]_ ^ Q(Ao + /i) - Q(Ao - /i) 

P 2sn(//) ^ ' 



Note that, for arbitrary number of interacting spins of the model, P and Q contain 
all the dynamical variables so that the Backlund maps touch all the spin sites. These 
maps associate to a given solution of the equations of motion fllip a new solution. By 
fixing the initial conditions, the generating function ([7]), and therefore the function 
7(A), is a constant independent of time. As will be clear in the next lines, if 7(A) is 
constant, then the Backlund transformations are actually rational maps (or, better to 
say, birational maps). The equation f lT6|) allow us to write the explicit transformations 
as follows: 



s 



+ 4 - PI - K)sl + + 4 - PI - al)sl - 2(a,/3, - ,,5,)sl) 
{zjal + 51- 4 - f5l)sl + + 4 + al + 5l)sl - 2.(/?,a, + ,,5u)sl) 
{{l^k<,k - Oik5k)sl + z(/3fc^-fc + ak6k)sl + {ak<.k + I^k5k)sl) 



(23) 



k 



where for brevity of notation we have introduced the functions (a^, 5^, A^, <jfc) de- 
fined by the following formulae: 

ttfc = sn(Afc - Ao) H cn(Afc - Aq) 

cn(^) 

(P + Q(2sn(/x)-PQ)) (P-g(2sn(^)-PQ)) ^^,^ 

= 2 + 2M^) ^^^^'^ ~ 

(P + Q(2sn(^)-FQ)) (P - g(2 sn(/i) - PQ)) ^ (24) 

= 2 2dii(;;i^ ^"^^'^ - 

q-fc = sn(Afc - Ao) -—^ cn(Afc - Aq) 

cn(/i) 

When the elliptic modulus k of the Jacobi elliptic functions is zero, the transformations 
(1231) coincide with those for the trigonometric Gaudin magnet as given in [12]. At this 
point we have to deal with the symplecticity of our maps. As the transformations are 
explicit, the direct path to prove their symplecticity might be to use ([6]) in order to 
show that indeed the Poisson structure is preserved. But this path is presumably not 
so plain because of the huge calculations. As in [T2j, we will follow a finer argument 
due to Sklyanin [2H]- Consider the relation f lT^ in an extended phase space, whose 
coordinate are given by (s^, s^, si, P, Q) and suppose that D{X) obeys to the quadratic 
Poisson bracket, as follows: 

{P'(A) ® 11, 11 ® D{t)} = [re(A - r), Z}(A) ® D{t)] (25) 
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In the extended space we have to re-define (fT2|) as: 



L{\)D{X- \o) = D{X- Xo)L{\) (26) 

In fact in the left hand side of the previous one has to use tilded variables also for -D(A) 
because (126!) defines the Backlund transformation in the extended phase space, where 
there is also a P and a Q. Note that in the new phase space the entries of D Poisson 
commutes with those of L. The key observation is that if both L and D have the same 
Poisson structure, given by equation (1251) . then this property holds true for LD and DL 
as well, because of the Poisson commutativity of the entries of L and D. This means 
that the transformation fl26l) defines a "canonical" transformation. Sklyanin showed 
[26] that if one now restricts the variables on the constraint manifold P = P and Q = Q 
the symplecticity is preserved; however this constraint leads to a dependence of P and 
Q on the entries of L, that for consistency must be the same as the one given by the 
equation ( 126|) on this constrained manifold. But there ( l26l) reduce to ( fT2i) . so that the 
map preserves the spectrum of L{X) and is canonical. What remains to show is that 
indeed (125!) is fulfilled by our -D(A). For the Backlund transformations of the rational 
Gaudin magnet the dressing matrix has the quadratic Poisson structure imposed by 
the rational r-matrix provided P and Q are canonically conjugated in the extended 
space [2S]- In the trigonometric case one need to have a non trivial bracket between 
P and Q in the extended space in order to guarantee the simplecticity of Backlund 
transformations [19]. As we will show in the next lines, in the elliptic case we found a 
non trivial bracket that, in the limit A; — )■ goes to the trigonometric result as given 
in [19]. By a direct inspection it is possible to show that ( 125|) entails the following 
brackets between the elements J^i, i = 0..3 [5]: 

{=-^1 '^o\ Jjk'S^i'^k (27) 



fc 



where (z, /c) is a cyclic permutation of (1, 2, 3) with J12 = /c^, ^23 = 1 — ^ J31 = — 1- 
With the following positions: 



(^0 



\2 



cn(//)dn(yu) 



sn(/i) (^1-2 sn{fi)PQ + P2g2 _ p (pg _ sn(/i))2 + cn(/i) 



cn(/i) 

tYi + ^2 dn(^) = PYo 
[ - ^2 dn(/i) = g(2 sn(^) - PQ) ^0 

(2^ 

after some calculations one can show that indeed (1271) are fuUfiUed provided that: 



f 1 - 2 sn(/i)PQ + P2Q2 _ p [(pg _ sn(/i))2 + cn(/i) 



2 (Q(PQ-2sn(M))-P)^ 
4 



cn(/i)dn(/i) 

(29) 
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So the symplecticity of our maps is proved. At this point let us to stress some remarks. 
The dressing matrix defined by the relations (128|) is, except for the multiplicative factor 
=5^0, completely equivalent to the dressing matrix as given by the equation (fT8|) . As 
explained before, given the homogeneity of the equation (I12p . a factor of proportionality 
between two dressing matrices is inessential as regards Backlund transformations, so by 
this point of view the definitions f l28p are compatible with f llSp . Secondly, the Poisson 
bracket fl29|) between P and Q reduces exactly, in the limit /c — )■ 0, to the bracket of 
the corresponding variables in the trigonometric case. 



5 Physical Backlund transformations 

In this section we will show that with an appropriate choice of the parameters Aq and 
/i in 0231) . the Backlund transformations map real solutions in real solutions, so in this 
sense the transformations can be considered "physical" . The choice amounts to require 
Aq to be a real number and /i to be a purely imaginary number. So, hereafter in this 
section, we put: 

fi = ie (Ao, e) G (30) 

The matrices Lj, j = 1..N defined in (fTTl) and corresponding to the real solutions 
(sj, s^,Sj) of the equations of motion are Hermitian. The request for physical Backlund 
transformations is equivalent to the request for Hermitian dressed matrices Lj. By (fT6|) 
we see that this means to have dressing matrices Dj proportional to unitary matrices. 
We claim that indeed when ( 130|) are fulfilled then Dj are of the form: 

where the bar means complex conjugation. For clarity let us make the following posi- 
tions: 

A+ = Ao + ie, A_ = A+ (32) 
We observe that, for the functions A, B, C, as defined in one has: 

A(A+) = A(A_), i?(A+) = C'(A„), C(A+) = S(A_). (33) 

These relations entails 7^(A+) = 7^(A_). Note that this last relation implies that 
the coefficients of the series of 7^ (A) with respect to A are real, consistently with the 
expansion ([7]). We recall that the matrices Dj are written in terms of P and Q, that 
are defined by the relations 

^ B{X^) A(A+) + 5(A+) Q(A+)-Q(A_) 

By specifying the sign of the function 7 on the Riemann surface by 7(A+) = — 7(A_), 
one has: 

- -oh 
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and this equation in turns implies that the matrices Dj = D{X) , with D{X) given 



by ( !T8|) . are of the form fl3T|) . with aj and (3j given by the following formulae: 

_ sn(z6)cn(A,)|Qp-l 
^ ^cn(^6)sn(A,)|g|2 + l 

sn(26) / Q + Q dn(Aj) Q - g ^ ^ ^ 

sn(A,) VigP + l ^ dn(2e) |gp + l. 

So, under the given assumptions, the matrices Dj are proportional to unitary matrices. 

6 Interpolating Hamiltonian flow 

Now we want to get the interpolating flow of the discrete dynamics generated by 
the maps 023]) • As we will see the Backlund transformation can be seen as a time 
discretization of a one-parameter (Aq) family of Hamiltonian flows with the difference 
2e playing the role of the time-step and with the Hamiltonian defining the interpolating 
flow given by 7(Ao), where 7(A) is defined in (120|) . 
First of all let we take the limit e — > 0. 
One has: 

Q= -"'°]-J''°' h-OM. (36) 

P = + O(t^). (37) 

7(Ao) 

One can carefully insert these expressions in the dressing matrix ( fTSl) to find: 

D{X - Ao) = 11 - 7 Do(A, Ao), (38) 

7(Ao)sn(A - Ao) 



where 

nn^^^^ ^(Ao)cn(A - Ao) B(Ao)±c(M + MAo)^dn(A - Ao 

^o(A, Ao) - |^ i.(Ao)+c(Ao) _ MAo)^dn(A - Ao) -A(Ao)cn(A - Ao) 

(39) 

In the limit e — the equation of the map LD = DL turns into the Lax equation for 
a continuous flow: 

L(A) = [L(A),M(A,Ao)]. (40) 
where the time derivative is defined as: 

L = hm (41) 

and the matrix M(A, Aq) is given by: 

M(A, Ao) = \ DoiX, Ao). (42) 

27(Ao)sn(A - Ao) 
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With the help of the Poisson brackets between the elements of the Lax matrix ([5]), the 
dynamical system (l40l) can be cast in Hamiltonian form: 

4.(A) = {H(Ao),L,,(A)}, e {1,2}, (43) 

with the Hamilton's function given by: 

•H(Ao) = 7(Ao) = VA^{Xo) + B{Xo)C{Xo). (44) 

So the Hamiltonian characterizing the interpolating flow is (the square root of) 
the generating function ([7j) of the whole set of conserved quantities. By choosing the 
parameter Aq to be equal to any of the poles (Aj) of the Lax matrix, the map leads to 
different maps {i?T'^*^}j=i..Ar, where BT^"^^ discretizes the flow corresponding to the 
Hamiltonian Hi, given by equation In fact, by posing Aq = 5 + Aj and taking the 
limit (5—7-0, the Hamilton's function (144|) gives: 

7(Ao) = y + - + 0(<5). (45) 
Si 

and the equations of motion take the form: 

4(A) = -{H„L,^{X)}, z,j G {1,2}. (46) 

Note that the corresponding interpolating Hamiltonian flows of the Backlund trans- 
formations for the trigonometric and rational Gaudin models can be obtained as the 
corresponding limiting case of this one. 



7 An application to the Clebsch model 

The Clebsch model [1] is an integrable case of the Kirchhoff equations [TT] describing 
the motion of a solid in an infinite incompressible fluid. 

If the solid has three perpendicular planes of symmetry and there are no external forces 
then the Kirchhoff system can be described in terms of an Hamiltonian (the kinetic 
energy of the system solid+fluid) quadratic and diagonal in the impulsive force p and 
impulsive pair J vectors, representing respectively the sum of the impulse and angular 
momentum of the solid and those applied by the solid to the boundary of the fluid in 
contact with it [T5] . 
The Hamiltonian then reads: 

T = l {a^ipn' + a,{pyf + a,{p^f) + \ {^.{rf + ^.{jyf + P,{rf) , (47) 

where Ui and Pi are a set of constants depending on the shape of the solid. Clebsch [1] 
discovered that if the following constraint on these quantities holds: 

^-—^ + + = 0, (48) 

P3 Pi P2 
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then the corresponding equations of motion are integrable. 

Since the dynamical variables are impulse and angular momentum vectors, the Lie- 
Poisson structure is defined by the e(3) algebra: 

{J\ J^} = e^^■'=J^ {J\p^} = e^^V, {P\P^} = 0. (49) 

where k belong to the set {x, y, z}. These brackets have two Casimirs: 

pj = ci, p^ = C2. (50) 

Now we will show how, by the means of a procedure of pole coalescence on the Lax 
matrix of the two-site elliptic Gaudin model, it is possible to obtain the Lax matrix 
for the Clebsch model governed by the following Hamiltonian: 

2n = {c- e{A + B)) {p^Y + c{pyf + (c + b{i - e)) {p^f + B{rf+ 

+ {Ak'' + B){jyf + {A + B){,rf. ^^^^ 

The construction of the integrability structure of the model is not new (cf. with 
[T7]). but here is briefiy reported for completeness. For more details about the pole 
coalescence procedure see [16], [17]. 

Note that the Clebsch's constraint f l48p holds for flHT]) . but here we have only four (not 
five) arbitrary constants, so actually we will obtain a special realization of f H7|) . Note 
also that in the case = one obtains the Kirchhoff top [21] . 

First let us introduce a contraction parameter, say e, and take in the Lax matrix 
([T]) Ai — 7- eAi and A2 — > eA2, where we recall that Ai and A2 are the two arbitrary 
parameters of the Gaudin model. By setting: 

J = si + S2, p = e(AiSi + A2S2) (52) 

and letting e — in ([T]) after this identification, one obtains the following expression: 

sn(A)L(A) = 

cn(A) r - idn(A) + ^^p^ - i^p^ 



sn(A) x^x^v^v;^ I gjj^;)^-) ^sn(A)J 

r + i dn(A) + ^-^^p^ + i^^py -cn(A) - 



(53) 

The overall term sn(A) multiplying L{X) on the l.h.s. of the previous equation can be 
obviously skipped, so in the following we assume that the Lax matrix of the model, 
that we will call L'^, is given only by the r.h.s. term of f l53|) . Note that by using (Q, it 
is readily seen that the variables J and p ( 152|) obey the Lie-Poisson algebra e(3) ( 149|) . 
The determinant of this matrix is the generating function of the integrals of motions. 
One has: 

- detmX)) = = H,- H,sn\X) + 2c, ^^^^^^^ + C2^ (54) 

sn(Aj sn'^^Aj 

where ci and C2 are the Casimirs fl50|) , while Hq and Hi are the two Poisson commuting 
integrals: 

Ho = {j^f + P {{jyf - {p^f) , Hi = 3' + {p^y - {{p^f + ip^Y) . (55) 

The physical hamiltonian f lSTjl is now obtained by the following linear combination: 

2n = AHo + BHi + Cci 
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7.1 Backlund transformations 



The ansatz for the dressing matrix for the model just described is inherited from that 
of the eUiptic Gaudin model since the r-matrix structure is preserved by the pole 
coalescence procedure. So for the general form of D{X) we again refer to ( fTSl) : 

pq — sn(/i) 

1 + iy ™(^) »^i+^2dn(A) 



D{X) 



3 sn{A) 
i.yi-,9'2dn{X) 
sn{A) 



sn(A) 

■J sn(A) 



with < 



cn(/i) 
z^i + ^2 dn(/i) = p 
z^i — ^2 dn(//) = g(2 sn{fi) — pq) 

(56) 

Exactly as for the Gaudin model, the variables p and q can be determined in terms of 
only one set of variables thanks to the spectrality property. The result is: 



q{\o + /^) 



L5,(A)-T(A) 



p 



g(Ao + /i) -q{\o- fi) 
2 sn(yu) 



(57) 



A=Ao+/i 

Again by choosing Aq real and fi purely imaginary one obtains transformations sending 
real variables into real variables. To write explicitly the maps one can, for example, take 
the residue at the pole in A = in the equivalence L'^{X)D{X — Aq) = D{X — Xo)L'^{X) 
and its value at A = K{k), where K{k) is the complete elliptic integral of the first kind 
fl^7|) . We write directly the real transformations by posing: 

^ = ie (Ao,e)eM^ 

They reads: 



(58) 



py 



r 



ft2 - 7i2 



2ho 

al + al + bl + bf 



p + 1- 



aobo + aobo ^ 
-P 







2h 



^py + i- 



4 + bl -^px _ . Qq^o - ao^o ^^ 



h 







-p- 



aobo + boao 



2ho 
p"" + i 



ho 



boao - aobo 



.py 



ho ho ho 

al + al- bj -]A jx ^ j^/ al-al-b\ + b\ _ ^, a^bi + aA 



(59) 



2hi 

al + al + b\ + bl 



2hi 



P 



jy ^ - Qi-Qi + ^i-^i jx _ . ai^i -ai&i ^. 



2hi 2 k' hi hi 

where the functions (cj, 6j, hi), i G {0, 1}, are defined by the following formulae: 

ao = sn(ze)cn(Ao) {\q\^ - l) - sn(Ao) (|g|' + l) 



bo = sn(ze) (g + g) - . dn(Ao) (g - g) 



ai 



cn(Ao 



dn{ie) 
|g|2+ 1) + A; 



, sn(ie)sn(Ao 



dn(Ao 

bl = sn(ze) (g + g) — k' 
ttiP + \b 



cn(ie)cn(Ao 
sn(ie) 



1) 



(60) 



hi 



dn(ie)dn(Ao) 

' ^e{o,l} 
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In the previous formulae k' is the complementary modulus of the Jacobi elliptic func- 
tions, k' = \fY--W\ the bar means complex conjugation. The expression for follows 
for example by the constraint J • p = J • p, but a cleaner expression can be found taking 
the value A = K{k) + \K'{k) in U{\)D{\ - \) = D{\ - Xo)L%X), where K'{k) is the 
complementary integral (l^7|) . It reads: 

jz ^ al-hl-al + bl _ ^ ^ ^3^2 + 02^3 ^1 ^^k ^H^!i±Mip?/ 
h h h 



where 



cn(ie)cn(Ao) 



03 = —ik sn(ie) in + q) 

• / (61) 
, sn le sn Ao , . ' ' 

dn(ie)cn(Ao) 



I, • 2 I 2 7,2 7,2 

/l = + ttg — O2 — O3 



Note that if in these transformations one poses A; = 0, then the Backlund transforma- 
tions for the Kirchhoff top as given in [21] are obtained. 

Now let us consider the interpolating Hamiltonian flow. We recall that the r-matrix 
structure of the model is that of Gaudin, so the Poisson brackets ([5]) works again by 
substituting A = Ll^, B = and C = L21. This is enough to ensure that the in- 
terpolating Hamiltonian is given, as for the Gaudin models, by the square root of the 
generating function of the integrals of the system (l54l) for A = Aq: 



uf\\ -rn ^ lu u 2f\\^o cn(Ao)dn(Ao) cn2(Ao) 

n{Xo) = T(Ao) = \lHi- iiosn^(Ao) + 2ci — h . (62) 

sni^AoJ sn {Aq} 

With a choice of the parameter Aq it is possible to obtain a discretization of the con- 
tinuous flow corresponding to each linear combination of the Hamiltonians Hq and Hi. 
A last remark: the symplecticity of these maps simply follows from the symplecticity 
of the maps for the ancestor Gaudin model, again thanks to the preservation of the 
r-matrix structure. 



8 Comments 

First of all we have to mention the lacks of our construction. Unlike the rational case 
(cf. with [H]), we are not able to give the generating function of the canonical transfor- 
mations defined by the maps f l23|) . However this isn't only a matter of technical difficul- 
ties; indeed for the xxx Gaudin model the two-parameters Backlund transformations 
can be written as the composition of two simpler one-parameter transformations: the 
same property holds true for the generating functions; yet in the trigonometric case 
a factorization of the dressing matrix cannot lead to a one parameter dressing matrix 
preserving all the symmetries of the problem ([20]). There are two possibilities: or 
one is able to find directly the two parameters generating function, or one should look 
for a symmetry-violating generating function such that their composition restores the 
symmetries. These details can be interesting by a quantum point of view for their 
potential connections with the Baxter's Q operator [77] . 
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In [21] it was shown how the discrete orbits defined by the Backlund transformations 
for the Kirchhoff top exactly interpolate, in some special cases and in general as a 
conjecture, the continuous orbits of the corresponding physical fiow, indicating that 
indeed the Backlund transformations can be considered an integrator (numerical or 
analytical) of the corresponding continuous differential problem. It could be interest- 
ing to understand how many of these results remain true for the transformations here 
given, both for Gaudin and Clebsch models. Works are in progress in this direction. 



A Notations and formulae 

Let Wi, W2 be complex numbers such that their ratio is not real and consider the lattice 
A generated by these numbers: 

A = {w E C : w = UiWi + n2W2, '^-2 ^ 

The WeierstraB zeta function is given by |2]: 

C(„, = iH.V('^ + i + 4) (63) 
u ^-^ \ u — W W I 



The WeierstraB p function is minus the derivative of C,: 

p{u) = -C{u) = ^ + E (r^-T^ - ^) (64) 

By denoting the period Wj, such that i^i + t(;2 + f^a = and defining the set e,, i = 1..3 
by ti = f?(^), then holds the relation [2]: 

p{uw) = 62 + ^4 M (^^) 

w^sn^-u, k) 

where tz7 = (ci — 62)"^ and the elliptic modulus for the Jacobi "sn" function is given by 
k'^ = ^^(63 — 62). The Jacobi elliptic functions sn(M, A;), cn(u, k) and dn(ti, k) satisfies 
the following quasi-periodic relations [2]: 

sn(u + 2mK + 2mA", A;) = (-l)™sn(M, k) 

cn{u + 2mK + 2inK\ k) = (-l)"^+"cn(M, A;) (66) 
dn{u + 2mK + 2inK' , k) = (-l)"dn(?/, A;) 

where K and K' are respectively the complete elliptic integral of the first kind and the 
complementary integral: 

K(k) = f '^^ = K\k) = f (67) 

Jo ^(l-t2)(l_pt2) v/(l-t2)(l_(l_A;2)t2) 

The following formulae are useful in proving ([7]): 

cn(x ± y) = cn(x)cn(y) =]= dn(a; ± t/)sn(x)sn(y) (68) 
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dn(x ±y) = dn(x)dn(?/) =]= k'^ciii(x ± y)sn(x)sn(y) (69) 

(70) 



w (C(rox) - C{ujy)) - ^^1^^ = a{y - x) 



a{x) = w (C(w(x)) - C(2ti7x)) - 

Equations ( l68l) and ( !69l) are only a consequence of addition formulae for the Jacobi ellip- 
tic functions, ( ITOl) can be proved in few lines. In fact suppose that x and y vary while y— 
X remains constant and equal to h. Differentiating /(x) = w {C,{wx) — C{w{x + h))) — 



sn(fe) 



sn(x)sn(x+6) 



with respect to x we see that this function is independent of x. In fact 

rn . 2/ , M ^ ( W sn(6)(sn(x)sn(x + 6))' 

f {x) = w [p[[x + h)w) - p[xw)) , , ^ , ;;72— 

(sn(x)sn(x + b)j 

By using the relation fl65|l and again the addition formulas for the Jacobi elliptic func- 
tions it is readily shown that f'{x) = 0, so f{x) is a constant, that we can take as 
a function of b. This implies the relation w {({wx) — ({zuy)) — s^"^^^^^-) = o,{y — x). 
By posing y = 2x in this equation we obtain the function a{x) as in fITOl) . Now let 
us consider closely the formula for —det{L{X)). For brevity we pose in the following 
Vi = X — Xi and Vij = Aj — Xj. From ([2]) we have: 



,j sn[Vi)sn[Vj) 



^ sn(wi)sn(^;j) 

„2 



(71) 



By adding and subtracting the quantities Xlj^j (dn(%)sf s| -|- A;^cn(wjj)sf sj) in the last 
equation and using ( l68l) and ( l69l) . we find: 



-rfet(L(A)) = 5^ _ ^ (5^.^^dn(z;,,) + k\\s]^n[v,,)) + 

sn I ) 

cn(t;ij)s,^s| sfsj + dvL{vij)s\s^j C^^) 



sn(f ^)sn(^^■) 



Now, using formula (170|) on the denominator of the last sum of equation (1721) and 
defining 



one reaches the result ([7]). 
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